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Abstract 

Magnetization process and adiabatic demagnetization of the antiferromagnetic Ising spin clusters with the shape of 
regular polyhedra (Platonic solids) are exactly examined within the framework of a simple graph-theoretical approach. 
While the Ising cube as the only unfrustrated (bipartite) spin cluster shows just one trivial plateau at zero magnetiza¬ 
tion, the other regular Ising polyhedra (tetrahedron, octahedron, icosahedron and dodecahedron) additionally display 
either one or two intermediate plateaux at fractional values of the saturation magnetization. The nature of highly de¬ 
generate ground states emergent at intermediate plateaux owing to a geometric frustration is clarified. It is evidenced 
that the regular Ising polyhedra exhibit a giant magnetocaloric effect in a vicinity of magnetization jumps, whereas the 
Ising octahedron and dodecahedron belong to the most prominent geometrically frustrated spin clusters that enable an 
efficient low-temperature refrigeration by the process of adiabatic demagnetization. 
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1. Introduction 

A recent progress in a targeted design of molecular 
nanomagnets, which afford paradigmatic examples of as¬ 
semblies of a finite number of interacting spin centers, 
has triggered a considerable interest to explore a mag¬ 
netic behavior of small spin clusters magnetically isolated 
from the environment BJi, L3.121, 14{] . Despite their simplic¬ 
ity, small spin clusters may still exhibit favorable mag¬ 
netic properties such as the quantum tunneling of mag¬ 
netization and slow spin relaxation observable in single¬ 
molecule magnets, which have an immense application 
potential in a development of new resources suitable for 
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quantum information processing and/or novel generation 
of high-density storage devices flfii. In addition, 
the molecular nanomagnets being composed of small spin 
clusters are ideal for an experimental testing of the limita¬ 
tions of physical theories especially when physical prop¬ 
erties of relevant model systems can be calculated without 
any uncontrolled approximation. 


Magnetic properties of small spin clusters are predom¬ 
inantly determined by the nature of magnetic interactions 
between the spin centers in a cluster, whereas the Heisen¬ 
berg superexchange coupling is usually the most domi¬ 
nant interaction term that basically influences characteris¬ 
tic features of molecular nanomagnets 0 , 000 ]. A lot 
of attention has been therefore paid to the study of antifer¬ 
romagnetic Heisenberg spin clusters, which may exhibit 
striking intermediate plateaux in low-temperature magne¬ 
tization curves that often macroscopically manifest high! 
non-trivial quantum ground states 
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14, 131HI 12 1 HI 19, 2(J • An influence of the exchange, 


dipolar and single-ion anisotropy on the low-temperature 
magnetization process of the antiferromagnetic Heisen¬ 
berg spin clusters has attracted much less attention so far 


112 1L 22, |23l 24]. For this reason, it appears worthwhile 


to investigate the low-temperature magnetization process 
of antiferromagnetic Ising spin clusters, which are also 
capable of displaying several intermediate magnetization 
plateaux on assumption that a spin cluster is geometrically 
frustrated Hi 2hll . 


The main goal of the present work is to examine the 
magnetization process and adiabatic demagnetization of 
the antiferromagnetic Ising spin clusters with the shape 
of regular polyhedra (Platonic solids), which surpris¬ 
ingly seem not to be dealt with previously. The zero- 
field thermodynamics of the regular Ising polyhedra with 
the uniform interaction have been explored in some de¬ 
tail by Syozi |27| and Fisher |28| within the framework 
of dual, decoration-iteration and star-triangle transforma¬ 
tions, while the regular Ising polyhedra with the mixed 
ferromagnetic and antiferromagnetic couplings have been 
examined in relation with the spin-glass physics 


31]. It is noteworthy that the competition between the an¬ 


tiferromagnetic order, spin frustration and magnetic held 
has been studied by Viitala and co-workers by considering 
antiferromagnetic Ising spin clusters with several clus¬ 
ter geometries as for instance octahedron, body centered 
icosahedron and cubooctahedron 1 32 , 33]. More recently, 
an exact enumeration of states has been employed in or¬ 
der to calculate numerically the residual entropy, magne¬ 
tocaloric properties and magnetization process of various 
geometrically frustrated Ising spin clusters composed of 
triangular units 00. 


The organization of this paper is follows. Exact re¬ 
sults for the partition function, free energy and magnetiza¬ 
tion of the regular Ising polyhedra are derived in Section 
[2] within the framework of a graph-theoretical approach. 
The most interesting results obtained for the magnetiza¬ 
tion process and adiabatic demagnetization of the regular 
Ising polyhedra are presented and discussed in detail in 
Section [3] Finally, some conclusions and future outlooks 
are drawn in Section [4] 


2. Regular Ising polyhedra 

Let us consider the Ising spin clusters with the ge¬ 
ometry of five regular polyhedra (tetrahedron, octahe¬ 
dron, cube, icosahedron and dodecahedron) as depicted in 
Fig -HI The regular Ising polyhedra can be defined through 
the following Hamiltonian 

N b N 

H = jY J SiSj-hY J S i , ( 1 ) 

< U > (=1 

where .S', = +1 represents the Ising spin placed at /th ver¬ 
tex of a regular polyhedron, the first summation accounts 
for the antiferromagnetic Ising-type exchange interaction 
J > 0 between adjacent spins, the second summation ac¬ 
counts for the Zeeman’s energy of individual magnetic 
moments in the external magnetic field h > 0 and finally, 
N (Nb) stands for the total number of vertices (edges) of 
a regular polyhedron that is simultaneously equal to the 
total number of spins (bonds). It is noteworthy that the 
considered Ising spin clusters are sufficiently small finite- 
size systems in order to calculate all physical quantities of 
interest quite rigorously. In the following, we will calcu¬ 
late the partition function, free energy and magnetization 
of the regular Ising polyhedra within a relatively simple 
graph-theoretical approach. 

The canonical partition function of the regular Ising 
polyhedra is defined through the relation 

Z = ^ exp(-/?#), (2) 

{5,1 

where /3 = 1 !(k K T), k K is Boltzmann’s constant, T is the 
absolute temperature and the summation 2(5 } is carried 
out over all possible configurations of a set of the Ising 
spins {Si}. While the interaction part of the Hamiltonian 
(the first term in Eq. 0])) depends on many specific details 
of a particular spin configuration, the Zeeman’s energy 
(the second term in Eq. (|TJ ) solely depends on the total 
spin St - E/Ii S , and is independent of any other specific 
detail of a given spin configuration. It is therefore quite 
convenient to sort available spin configurations accord¬ 
ing to the total spin S y, because only the lowest-energy 
spin configuration from a given set of microstates with the 
same total spin S t may eventually become a ground state. 
In addition, it is sufficient to consider only spin configu¬ 
rations with non-negative values of the total spin S t > 0, 
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Figure 1: The Ising spin clusters with the shape of five regular polyhedra. The upper panel displays a steric arrangement of the regular polyhedra, 
while the lower panel shows the equivalent planar lattice graphs. 


since the exchange energy is invariant under the reversal 
of all Ising spins (5,- —> -.S', for each i) and the Zeeman’s 
energy merely changes its sign. 

To proceed further with a calculation, let us exploit the 
one-to-one correspondence between the Ising spins S ; = 
+ 1 and the newly defined two-valued variables n,- = 0,1 
established through the formula 

St = l-2m, (3) 


which connects the ’up’ spin state Si = 1 to the value 
n, = 0 and the ’down’ spin state Sj - -1 to the other 
value n, = 1 . The Hamiltonian (Q} can be consequently 
rewritten in terms of newly defined two-valued variables 
n, to the following form 


H = J 


N b -2j] (rn +rij- 2n,n ^ 

<U> 

N 

V-2 


(4) 


;=i 


It is quite evident from Eq. Q that the exchange energy 
N b J of the fully aligned spin configuration can be lowered 
by an amount -27 per each unequally oriented pair of ad¬ 
jacent Ising spins, but each Ising spin flipped in opposite 
to the external magnetic field 5, = -1 (n, = 1) simultane¬ 
ously raises the Zeeman’s energy —Nh by an amount 2h. 
Apparently, the lowest-energy spin configuration (ground 


state) at a given magnetic field must be a reasonable com¬ 
promise that would minimize a sum of both parts of the 
overall energy. 

The alternative form of the Hamiltonian (QJ also sug¬ 
gest an interesting mapping correspondence between the 
summation over all available configurations of the Ising 
spins and a relatively simple graph-theoretical counting 
problem. When the Ising spins oriented in opposite to 
the magnetic field S, = -1 («,■ = 1) are identified as 
vertices of the equivalent planar graph (see lower panel 
in Fig. |T}, the combinatorial problem of finding all possi¬ 
ble spin configurations becomes equivalent to the problem 
of finding all induced subgraphs of a given planar graph. 
As a matter of fact, the total number of the flipped Ising 
spins (S, = -1 , rij = 1) then directly equals to the total 
number of vertices n t = XZi n ’ °f the corresponding in¬ 
duced subgraph and, respectively, the total number of un¬ 
equally oriented adjacent spin pairs ( n i + n j ~ 2m n j) 
equals to a sum over complementary degrees of all ver¬ 
tices present in a given induced subgraph d, = E/Ii = 
YjZi n i(d - di) (d is degree of a vertex in a full graph and 
d, is degree of the /th vertex in a given induced subgraph). 
Hence, it follows that the overall energy of each individ¬ 
ual spin configuration can be straightforwardly calculated 
from the total number of vertices n, and the sum of all 
complementary degrees d t of the corresponding induced 
subgraph according to the formula 

E = J(N h - 2d t ) - h(N - 2n t ). (5) 
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Table 1: Spin configurations of the Ising tetrahedron classified accord¬ 
ing to the total spin St > 0, the total number of flipped Ising spins 
(ill), the degeneracy (deg), the total number of unequally oriented ad¬ 
jacent spin pairs (d t ), the overall energy and the corresponding induced 
subgraph (see Fig.[2]for schematic representation of induced subgraphs). 
The values n, and d t coincide with the total number of vertices in a given 
subgraph and the sum of their complementary degrees, respectively. 


S T 

n, 

deg 

d t 

energy 

subgraph 

4 

0 

1 

0 

6 J - 4 h 

0 

2 

1 

4 

3 

—2h 

1A 

0 

2 

6 

4 

-2 J 

2B 


It is worthwhile to remark that the graphs corresponding 
to the regular Ising polyhedra are vertex-transitive, which 
substantially simplifies the relevant calculation of com¬ 
plementary degrees of vertices in respective induced sub¬ 
graphs. 

Let us illustrate the aforedescribed calculation proce¬ 
dure on particular case of the regular Ising polyhedra 
starting from the simplest example, which is the Ising 
tetrahedron (see Tab. H}. When all four spins are in the 
’up’ spin state .S', = 1 (n, = 0) one obtains a fully aligned 
spin configuration of the Ising tetrahedron with the total 
spin St - 4, which corresponds according to the estab¬ 
lished mapping equivalence to a very special subgraph 0 
(empty graph) with a zero vertex set. All other microstates 
of the Ising tetrahedron with St > 0 can be then obtained 
from this fully aligned spin configuration by turning a de¬ 
sired number of the Ising spins to their ’down’ spin state, 
whereas it is sufficient to flip no more than a half of the 
total number of spins owing to the invariance of the ex¬ 
change energy with respect to the inter-change .S', —> -S 
Bearing this in mind, only two additional spin configura¬ 
tions of the Ising tetrahedron should be taken into account 
with either one or two inverted spins, whereas the former 
spin configuration with St = 2 corresponds to an induced 
subgraph containing an isolated vertex (the subgraph 1A 
in Fig. [3 and the latter spin configuration with St - 0 
to an induced subgraph involving two vertices incident to 
a common edge (the subgraph 2B in Fig. [3- It is worth¬ 
while to remember that the overall energy of spin config¬ 
urations with St < 0 can be obtained from that ones with 
S t > 0 simply by inverting the sign of Zeeman’s term 
(the exchange energy is invariant under this reversal). 

The aforedescribed calculation procedure can be 
straightforwardly adapted also to the larger Ising spin 
clusters, so let us merely emphasize the most essential dif¬ 


Table 2: The Ising octahedron (see Tab[l]for notation). 


S T 

n, 

deg 

d t 

energy 

subgraph 

6 

0 

1 

0 

12J - 6h 

0 

4 

1 

6 

4 

4 J - 4 h 

1A 

2 

2 

3 

8 

1 

4 - 

1 

2A 

2 

2 

12 

6 

-2/2 

2B 

0 

3 

12 

8 

-4/ 

3C 

0 

3 

8 

6 

0 

3D 


ferences to emerge with increasing a spin-cluster size. It 
can be already understood from the spin configurations of 
the Ising octahedron listed in Tab.|3that there may appear 
a few different spin configurations with the same value of 
the total spin S but which need not have the same con¬ 
figurational energy. For instance, two overturned spins 
in the Ising octahedron may be or not nearest neighbours 
with respect to each other and hence, the relevant spin 
configurations necessarily correspond to two different in¬ 
duced subgraphs: one is being composed of two vertices 
incident to a common edge (the subgraph 2B in Fig. [3 
and the other one consists of two isolated vertices (the 
subgraph 2A in Fig. [3- It is quite clear that the latter spin 
configuration has lower energy due to a greater total num¬ 
ber of unequally oriented adjacent spin pairs (i.e. greater 
sum of complementary degrees of vertices in a given in¬ 
duced subgraph) in comparison with the former spin con¬ 
figuration, which can be consequently excluded from pos¬ 
sible candidates for a ground state. A similar finding con¬ 
cerns with two spin configurations of the Ising octahedron 
with three reversed spins, which may be or not all nearest 
neighbours with respect to each other (see the subgraphs 
3C and 3D in Fig. [3. For easy orientation, the spin con¬ 
figurations that represent actual ground state over a finite 
range of magnetic fields are highlighted bold in Tab.[l]||] 
All spin configurations of the Ising cube with the non¬ 
negative total spin St >0 are listed in Tab. [3 As one can 
see, a greater capability of the Ising cube in generating 
the spin configurations with the same total spin S t can be 
responsible for an accidental degeneracy of a few differ¬ 
ent microstates with the identical configurational energy. 
Another interesting observation is that the lowest-energy 
spin configurations from a given set of microstates with 
the same value of the total spin S t are the ones, in which 
occupation of the nearest-neighbour vertices by any pair 
of inverted spins is forbidden because the overall energy 
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Figure 2: A schematic representation of the induced subgraphs, which correspond to all possible spin configurations of the Ising tetrahedron, 
octahedron, cube and icosahedron, as well as, the lowest-energy spin configurations of the Ising dodecahedron. 
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Table 3: The Ising cube (see Tab[I]for notation). 


S T 

n, 

deg 

d t 

energy 

subgraph 

8 

0 

1 

0 

1 

00 

0 

6 

1 

8 

3 

6 / - 6/2 

1A 

4 

2 

16 

6 

-4/2 

2A 

4 

2 

12 

4 

4/ - 4/2 

2B 

2 

3 

8 

9 

-s: 

<N 

1 

so 

1 

3A 

2 

3 

24 

7 

-s: 

<N 

l 

<N 

1 

3B 

2 

3 

24 

5 

2J - 2/2 

3C 

0 

4 

2 

12 

-12/ 

4A 

0 

4 

6 

8 

-4 J 

4C 

0 

4 

24 

8 

-4 J 

4D 

0 

4 

24 

6 

0 

4E 

0 

4 

8 

6 

0 

4G 

0 

4 

6 

4 

4/ 

4H 


is always minimized by the induced subgraphs entirely 
composed of a set of isolated vertices (see the subgraphs 
1 A, 2A, 3A and 4A in Fig. 0. This special property is a 
direct consequence of the fact that the Ising cube (as the 
only regular polyhedron) is a bipartite spin cluster. 

Let us also shortly comment on a striking accidental de¬ 
generacy of the lowest-energy spin configurations of the 
Ising icosahedron, which constitute the highly degener¬ 
ate ground-state manifold with the total spin St - 0 (see 
Tab. 0. It actually turns out that the true ground state 
of the Ising icosahedron at small enough magnetic fields 
is the highly degenerate ground-state manifold, in which 
six overturned spins either form a path of six vertices (the 
subgraph 6B in Fig. [2J or one isolated vertex plus a cycle 
with five vertices (the subgraph 6C in Fig. 0. The other 
highly degenerate lowest-energy spin configurations of 
the Ising icosahedron can be found among the microstates 
with five inverted spins (i.e. the total spin St - 2), which 
either constitute two disjoint paths composed of two and 
three vertices (the subgraph 5B in Fig. 0 or an isolated 
vertex supplemented with a path graph of four vertices 
(the subgraph 5C in Fig. 0. It is worthwhile to remark, 
however, that the lowest-energy spin configurations with 
the total spin St = 2 do not represent true ground states 
unlike the spin configurations with the total spin St = 0. 

It is quite evident that the procedure illustrated in above 
on relatively small regular Ising spin clusters with the to¬ 
tal number of spins N < 12 would become rather cum¬ 
bersome when applying it to much larger spin clusters as 


Table 4: The Ising icosahedron (see Tab[I|for notation). 


S T 

n, 

deg 

d, 

energy 

subgraph 

12 

0 

1 

0 

30/ - 12/i 

0 

10 

1 

12 

5 

20 J - 10/; 

1A 

8 

2 

36 

10 

10/ - 8/; 

2A 

8 

2 

30 

8 

14/ - 8 h 

2B 

6 

3 

20 

15 

—6h 

3A 

6 

3 

120 

13 

4/ - 6/; 

3B 

6 

3 

60 

11 

so 

1 

oo 

3C 

6 

3 

20 

9 

12/ - 6 h 

3D 

4 

4 

30 

18 

-6 J - 4 h 

4B 

4 

4 

75 

16 

-s; 

1 

<N 

1 

4C 

4 

4 

120 

16 

-s; 

1 

<N 

1 

4D 

4 

4 

120 

14 

2/ - 4 h 

4E 

4 

4 

60 

14 

2/ - 4/i 

4F 

4 

4 

60 

12 

6/ - 4/7 

41 

4 

4 

30 

10 

10/ - 4/7 

4J 

2 

5 

60 

19 

-s; 

<N 

1 

OO 

1 

5B 

2 

5 

60 

19 

-s; 

<N 

1 

OO 

1 

5C 

2 

5 

180 

17 

-s; 

<N 

1 

1 

5D 

2 

5 

60 

17 

-s; 

<N 

1 

1 

5E 

2 

5 

60 

17 

-s; 

(N 

1 

1 

5F 

2 

5 

60 

15 

-2/z 

5G 

2 

5 

12 

15 

-2/z 

5H 

2 

5 

120 

15 

-2 h 

51 

2 

5 

60 

15 

-2h 

5J 

2 

5 

60 

13 

4/ - 2h 

5K 

2 

5 

60 

11 

1 

OO 

5L 

0 

6 

60 

20 

-10/ 

6B 

0 

6 

12 

20 

-10/ 

6C 

0 

6 

40 

18 

-6/ 

6D 

0 

6 

10 

18 

-6/ 

6E 

0 

6 

30 

18 

-6/ 

6F 

0 

6 

20 

18 

-6/ 

6G 

0 

6 

120 

18 

-6/ 

6H 

0 

6 

120 

18 

-6/ 

61 

0 

6 

60 

16 

-2/ 

6J 

0 

6 

60 

16 

-2/ 

6K 

0 

6 

30 

16 

-2/ 

6L 

0 

6 

30 

16 

-2/ 

6M 

0 

6 

120 

16 

-2/ 

6N 

0 

6 

120 

14 

2/ 

60 

0 

6 

60 

12 

6/ 

6P 

0 

6 

20 

12 

6/ 

6Q 

0 

6 

12 

10 

10/ 

6R 
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Table 5: The Ising dodecahedron (see Tab[I]for notation). Only the 
lowest-energy spin configurations with 5 r > 0 are listed. 


S T 

n, 

deg 

d, 

energy 

subgraph 

20 

0 

1 

0 

307 - 20/z 

0 

18 

1 

20 

3 

247- 18/z 

1A 

16 

2 

160 

6 

187- 16/z 

2A 

14 

3 

660 

9 

127 - 14/z 

3A 

12 

4 

1510 

12 

67 - 12h 

4A 

10 

5 

1912 

15 

- 10 /z 

5A 

8 

6 

1240 

18 

1 

Os 

1 

OO 

6 A 

6 

7 

320 

21 

-127 - 6 /z 

7A 

4 

8 

5 

24 

1 

h- 1 
00 

1 

4- 

*5- 

8A 

2 

9 

840 

23 

-167 - 2/z 

9A 

2 

9 

60 

23 

-167 - 2/z 

9B 

0 

10 

240 

24 

-187 

10A 


for instance the Ising dodecahedron with N = 20. It is 
should be nevertheless pointed out that the main difficulty 
by adopting this calculation procedure for more exten¬ 
sive Ising spin clusters does not concern with enumerat¬ 
ing types of induced subgraphs inherent to individual spin 
configurations (this task is relatively straightforward), but 
rather in a combinatorial assignment of respective degen¬ 
eracies that must be attributed to each individual induced 
subgraph. However, the mathematical procedure worked 
out previously purely on analytic grounds is amenable 
for the numerical implementation as well. To illustrate 
the case, we have implemented a set of procedures in 
MAPLE computer algebra system in order to enumerate 
all induced subgraphs (and their respective occurrences), 
which are pertinent to available spin configurations of the 
Ising dodecahedron. For simplicity. Tab. [5] summarizes 
just the lowest-energy spin configurations of the Ising 
dodecahedron for non-negative values of the total spin 
Sr > 0. It can be readily understood from Tab.[5]that the 
spin configuration with eight disjunct flipped spins (the 
induced subgraph 8 A in Fig. 0 is favorable ground state 
of the Ising dodecahedron in a rather wide interval of the 
magnetic fields ranging from zero up to saturation field. 

At this stage, the overall energy spectrum of the reg¬ 
ular Ising polyhedra obtained within the simple graph- 
theoretical approach can be substituted to the definition 
of the partition function (0 in order to get full thermody¬ 
namics. For completeness, let us merely quote the final 
exact expressions for the partition function of the Ising 


tetrahedron 


Zj = 6 exp (2 j3J) + 8 cosh(2y3/z) 

+ 2exp(-6//7)cosh(4///z), ( 6 ) 

the Ising octahedron 

Z 0 = 8 + 12 exp(4 J3J) 

+ 6cosh(2///z)[exp(4/37)+4] 

+ 12exp(-4/)7)cosh(4///z) 

+ 2exp(-12//7)cosh(6/l/z), (7) 

the Ising cube 


Z c = 2 exp(12/)7) + 30 exp(4//7) + 32 + 6 exp(-4/37) 

+ 16 cosh(2y8/z)[exp (6/3 J) + 6 cosh (2/3 J)] 

+ 8 cosh(4/)/z)[4 + 3 exp(-4/37)] 

+ 16 exp (- 6 /k/)cosh( 6/5/z) 

+ 2 exp(-12//7) cosh( 8 /?/z), ( 8 ) 


and the Ising icosahedron 

Z/ = 72 exp (10/17) + 340 exp ( 6/3J ) + 300 exp (2/3J) 

+ 120 exp (-2/3J) + 80 exp (-6/3J) + 12 exp (-10 f3J) 
+ 24 cosh(2/l/z) [ 10 exp( 8/3 J) + 25 exp(4 /3J) + 21] 

+ 120cosh(2/)/z)[exp(-4/17) + exp(-8/I7)] 

+ 30 cosh(4/f/z)[2 exp (6/3 J) + 13 exp (2/3 J)] 

+ 120 cosh(4/)/z)[3 exp (-2/3 J) + exp ( -6/3J )] 

+ 60 exp (-10/17) cosh(4/?/z) 

+ 40exp(-12/37)cosh(6/l/z) 

+ 40 cosh(6/l/z)[l + 6 exp(-4 f3J) + 3 exp(-8/17)] 

+ 12 cosh(8/l/z)[6 exp (-10/17) + 5 exp (-14/17)] 

+ 24 exp( -20/17) cosh( 10/1/z) 

+ 2exp(-30/17)cosh(12/l/z). (9) 


Gibbs free energy normalized per one spin can be directly 
calculated from the partition functions of the regular Ising 
polyhedra according to the relation 

g = -^k B T\nZ, (10) 


which subsequently allows a straightforward calculation 
of the magnetization per spin 


dg 1 <9 In Z 
dh N d(J3h) 


(ID 
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as well as the entropy per spin 


tetrahedron 


s = _5* = ^(,„ z-i^V 

dT N\ p dp ) 


( 12 ) 


For brevity, let us merely quote the final exact analyti¬ 
cal formula for the single-site magnetization of the Ising 
tetrahedron 

2 

m T = — [2 sinh(2y8/z) + exp (-6 pj) sinh(4 ph )], (13) 

Z T 


the Ising octahedron 
2 

m 0 = — {sinh(2y6/z)[exp(4y67) + 4] (14) 

Zo 

+ 4 exp(-4y67) sinh(4y8/z) + exp(-12y67) sinh(6y8/z)}, 
the Ising cube 
2 

me = —{2 sinh(2yS/0[exp(6/l/) + 6 cosh(2y67)] 

Z c 

+ 2sinh(4j8/r)[3exp(-4/37)+4] (15) 

+ 6 exp (-6 pj) sinh(6 ph) + exp(-12 pj) sinh(8 ph )}, 

and the Ising icosahedron 

2 

mi = — |2 sinh(2/l/z)[10exp(8/l./) + 25 exp(4 pj) + 21] 
Z/ 

+ 10 sinh(2/I/!)[exp(-4y6/) + exp(-8yS7)] 

+ 5 sinh(4/I/z)[2 exp (6 pj) + 13 exp (2 pj)] 

+ 20 sinh(4/l/z)[3 exp (-2 pj) + exp (-6 pj)] 

+ 10exp(-10y67)sinh(4/l/2) 

+ 10 exp( -12 pj) sinh(6 ph) 

+ 10 sinh(6y6/i)[l + 6 exp(-4 pj) + 3 exp(-8yS7)] 

+ 4sinh(8y6/z)[6exp(-10/l/) + 5 exp (-14 pj)] 

+ 10 exp(-20/6/) sinh(10yS/z) 

+ exp(-30j6/)sinh(12/I/z)}. (16) 

It should be noted that the final expressions for the en¬ 
tropy of the regular Ising polyhedra obtained according to 
Eq. GU) are too cumbersome to write them down here 
explicitly. 
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Figure 3: The upper panel: 3D surface plot for the magnetization of the 
Ising tetrahedron against the magnetic field and temperature; the lower 
panel: adiabatic temperature changes of the Ising tetrahedron upon vary¬ 
ing the magnetic field. The magnetization is normalized with respect to 
its saturation value and entropy per spin. 



3. Results and discussion 

In this part, let us explore in detail the magnetization 
process and magnetocaloric properties of the regular Ising 
polyhedra. In particular, we will investigate a microscopic 
nature of highly degenerate ground states, which are man¬ 
ifested in low-temperature curves as profound intermedi¬ 
ate plateaux at fractional values of the saturation magne¬ 
tization. 

The magnetization of the Ising tetrahedron is plotted in 
the upper panel of Fig. [3] against temperature and mag¬ 
netic field. It can be clearly seen from this figure that a 
zero-temperature magnetization curve of the Ising tetra- 
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hedron shows two intermediate plateaux at zero and one- 
half of the saturation magnetization. The magnetization 
is zero in a field range h/J e (0,1), which favors the six¬ 
fold degenerate ground state with two spins oriented in 
a direction of the magnetic field and the other two spins 
pointing in opposite direction. One out of two spins ori¬ 
ented contrary to the magnetic field is flipped at the first 
critical field hj J = 1, above which the four-fold degen¬ 
erate ground state emerges with three spins aligned into 
the magnetic field and one spin pointing in opposite di¬ 
rection. The aforementioned ground state manifest it¬ 
self as an intermediate plateau at one-half of the satura¬ 
tion magnetization, which can be detected in a field range 
h/J e (1,3). Finally, the fully polarized state with all 
four spins aligned towards the magnetic field becomes the 
ground state above the second critical field h c /J = 3. It 
should be pointed out that exact magnetization jumps and 
plateaux occur only at zero temperature, because a ris¬ 
ing temperature generally causes a gradual smoothening 
of the magnetization curve. As a result, magnetization 
plateaux generally shrink and magnetization jumps are 
smeared out with increasing temperature. 

Let us also examine magnetocaloric properties of the 
Ising tetrahedron, which can be particularly interesting 
especially in a vicinity of the aforedescribed magnetiza¬ 
tion jumps. For this purpose, the lower panel of Fig. [3] 
displays a few isentropy lines of the Ising tetrahedron in 
the magnetic field-temperature plane, which can be alter¬ 
natively viewed as a temperature response with respect 
to a variation of the magnetic field under the adiabatic 
condition. It actually turns out that the Ising tetrahedron 
may exhibit an enhanced magnetocaloric effect in a vicin¬ 
ity of two magnetization jumps, above (below) which an 
abrupt rise (drop) in temperature is achieved by decreas¬ 
ing magnetic field whenever the entropy per spin is set 
close enough to the ground-state degeneracy pertinent to 
the first and second critical held s/k b = |lnl0 = 0.576 
and S/fa = | In 5 = 0.402, respectively. 

In contrast to the Ising tetrahedron, the magnetization 
curve of the Ising octahedron shows just one wide plateau 
at one-third of the saturation magnetization in a respective 
held range h/J e (0,4), see the upper panel in Fig. [4] It is 
worthwhile to remark that the intermediate plateau at zero 
magnetization is totally absent in a relevant magnetization 
process, because three spin configurations with two an- 
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Figure 4: The upper panel: 3D surface plot for the magnetization of the 
Ising octahedron against the magnetic field and temperature; the lower 
panel: adiabatic temperature changes of the Ising octahedron upon vary¬ 
ing the magnetic field. The magnetization is normalized with respect to 
its saturation value and entropy per spin. 


tipodal spins oriented against the magnetic held are ener¬ 
getically favorable before any spin conhguration with the 
total spin S j — 0 as soon as the magnetic held becomes 
non-zero. In virtue of this fact, the zero-temperature mag¬ 
netization curve of the Ising octahedron exhibits an abrupt 
jump from zero to one-third of the saturation magnetiza¬ 
tion right at zero held, whereas the intermediate one-third 
plateau terminates at the saturation held h c /J = 4 that de¬ 
termines a held-induced transition towards the fully sat¬ 
urated state. The effect of temperature upon a gradual 
smoothening of the magnetization curve is quite similar 
as in the case of the Ising tetrahedron. 

Adiabatic changes of temperature achieved upon vary- 
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ing the magnetic field are depicted for the Ising octa¬ 
hedron in the lower panel of Fig. [4] Apparently, the 
Ising octahedron shows an enhanced magnetocaloric ef¬ 
fect near the saturation field provided that the entropy 
per spin is close to the ground-state degeneracy relevant 
for the saturation field s/k b = g In 10 = 0.384. How¬ 
ever, it is worth noticing that an absence of interme¬ 
diate plateau at zero magnetization has significant im¬ 
pact on magnetocaloric properties of the Ising octahe¬ 
dron. The absence of zero magnetization plateau bears 
a close relation to a giant magnetocaloric effect observ¬ 
able near zero magnetic field. If the entropy is selected 
close enough to the ground-state degeneracy emergent at 
zero field s/k B = i In 18 = 0.482, then, the adiabatic de¬ 
magnetization of the Ising octahedron is characterized by 
a rapid decrease in temperature as the magnetic field di¬ 
minishes. It should be stressed that the isentropy line at 
j/^b = | In 18 = 0.482 reaches the absolute zero temper¬ 
ature with an infinite gradient as the magnetic field tends 
to zero, which suggest the fastest possible cooling process 
that would be quite superior with respect to that of param¬ 
agnetic salts. A refrigeration with the help of adiabatic de¬ 
magnetization of the Ising octahedron might be therefore 
of substantial technological importance, since the efficient 
cooling persists in a rather wide range of the magnetic 
fields and the ultra-low temperatures can be consequently 
reached quite efficiently. 

The magnetization curve of the Ising cube as the only 
non-frustrated regular polyhedron consists of a single 
plateau at zero magnetization, which is subsequently fol¬ 
lowed by a direct magnetization jump from zero to the 
saturation magnetization (see the upper panel in Fig. 0. 
The bipartite nature of the Ising cube enables an existence 
of two-fold degenerate antiferromagnetic ground state at 
sufficiently low magnetic fields with four reversed spins, 
which undergo a simultaneous spin flip to the fully po¬ 
larized state at the saturation field h c /J = 3. It is note¬ 
worthy that the magnetization curve becomes smoother 
upon increasing temperature and there is a huge degener¬ 
acy precisely at the saturation field, at which all lowest- 
energy spin configurations with S j > 0 coexist together. 
It could be concluded that the vertex-transitive bipartite 
Ising spin clusters cannot exhibit any other magnetization 
plateau except that at zero magnetization. 

The adiabatic demagnetization of the Ising cube is dis¬ 
played in the lower panel of Fig. 0 which shows a con- 
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Figure 5: The upper panel: 3D surface plot for the magnetization of the 
Ising cube against the magnetic field and temperature; the lower panel: 
adiabatic temperature changes of the Ising cube upon varying the mag¬ 
netic field. The magnetization is normalized with respect to its saturation 
value and entropy per spin. 

tour plot of entropy as a function of temperature and mag¬ 
netic field. The giant magnetocaloric effect is observed 
in a vicinity of the saturation field on account of a huge 
degeneracy of the lowest-energy spin configurations with 
St > 0 on assumption that the entropy per spin is set to 
s/ks = yr In 35 = 0.296. However, a rather steep de¬ 
crease in temperature observable slightly above the satu¬ 
ration field is successively followed by a subsequent steep 
increase in temperature below the saturation field, which 
consequently precludes a possible utilization of a giant 
magnetocaloric effect for cooling purposes. 

The most remarkable magnetization process evidently 
exhibits the Ising icosahedron, which displays in total 
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Figure 6: The upper panel: 3D surface plot for the magnetization of 
the Ising icosahedron against the magnetic field and temperature; the 
lower panel: adiabatic temperature changes of the Ising icosahedron 
upon varying the magnetic field. The magnetization is normalized with 
respect to its saturation value and entropy per spin. 


three different intermediate plateaux at zero, one-third 
and one-half of the saturation magnetization (see the up¬ 
per panel in Fig. [6j- At sufficiently low magnetic fields 
h/J e (0,1) the ground state is highly (72-fold) degener¬ 
ate due to an accidental degeneracy of two lowest-energy 
spin configurations, in which six reversed spins are lo¬ 
cated at vertices that either form an open chain of six 
sites or a closed chain of five sites plus one additional iso¬ 
lated site (see the corresponding induced subgraphs 6B 
and 6C in Fig. [2}. The intermediate one-third magneti¬ 
zation plateau develops at moderate fields h/J e (1,3) on 
account of the another highly (30-fold) degenerate ground 
state, in which two out of four reversed spins are located 


at isolated vertices and the other two inverted spins at ad¬ 
jacent vertices (see the corresponding induced subgraph 
4B in Fig. |2}. The last intermediate plateau at one-half 
of the saturation magnetization occurs in the field range 
li/Je (3,5), which prefers the highly (20-fold) degener¬ 
ate ground state with three isolated spins oriented in oppo¬ 
site to the magnetic field (see the corresponding induced 
subgraph 3A in Fig. [2}. Apparently, the high degeneracy 
of all aforedescribed ground states is closely connected 
to a geometric spin frustration of the Ising icosahedron, 
which is progressively lowered with an increase of the ex¬ 
ternal magnetic held until the degeneracy is completely 
lifted above the saturation held hJJ = 5 due to the unique 
fully ordered spin conhguration. 

To achieve an enhanced magnetocaloric effect during 
the adiabatic demagnetization of the Ising icosahedron 
(the lower panel in Fig. ©, it is necessary to set entropy 
close enough to an appropriate ground-state degeneracy 
at a critical held related to the respective magnetization 
jump. Thus, a gradual decrease in the magnetic held is 
responsible for the fast cooling (heating) slightly above 
(below) the relevant critical held only if entropy per spin 
is close to s/k b = pr In 222 = 0.450 nearby the hrst crit¬ 
ical held hJJ = 1, s/kn = In 50 = 0.326 nearby the 
second critical held hJJ = 3 and s/kn = p In 69 = 0.353 
nearby the saturation held hjJ = 5. 

The magnetization process of the Ising dodecahedron 
(the upper panel in Fig. |7j is quite reminiscent of that 
of the Ising octahedron. The zero-temperature magneti¬ 
zation curve of the Ising dodecahedron shows an abrupt 
jump at zero magnetic held, which is consecutively fol¬ 
lowed by a rather wide intermediate plateau at one-hfth 
of the saturation magnetization ending up only at the sat¬ 
uration held hJJ = 3. The ground state of the Ising do¬ 
decahedron is just hve-fold degenerate despite of a strong 
spin frustration in a low-held region h/J e (0,3), where 
hve different spin conhgurations with eight isolated spins 
aligned contrary to the magnetic held become the lowest- 
energy microstates. The magnetization plateau and jump 
are gradually smoothened with an increase in temperature 
until they completely vanish above a certain temperature. 

One of the most important features of the magneti¬ 
zation process of the Ising dodecahedron is the absence 
of zero magnetization plateau, which results from an in¬ 
capability of microstates with the total spin St = 0 
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4. Conclusion 
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Figure 7: The upper panel: 3D surface plot for the magnetization of 
the Ising dodecahedron against the magnetic field and temperature; the 
lower panel: adiabatic temperature changes of the Ising dodecahedron 
upon varying the magnetic field. The magnetization is normalized with 
respect to its saturation value and entropy per spin. 


or 1 to further reduce the configurational energy of the 
ground-state manifold composed of five microstates with 
the total spin S j = 2. Owing to this fact, the giant 
magnetocaloric effect can be detected during the adia¬ 
batic demagnetization of the Ising dodecahedron in a very 
low-field region on assumption that entropy per spin is 
close enough to the ground-state degeneracy at zero field 
s/ks = ^jln250 = 0.276. From this perspective, the 
Ising dodecahedron could be regarded together with the 
Ising octahedron as the most prominent frustrated Ising 
spin cluster that would allow efficient cooling through the 
adiabatic demagnetization. 


The present work deals with the magnetization process 
and magnetocaloric properties of the regular Ising polyhe- 
dra (tetrahedron, octahedron, cube, icosahedron and do¬ 
decahedron), which are rigorously calculated by means 
of a relatively simple graph-theoretical approach. Apart 
from the exact analytic calculations, the proposed proce¬ 
dure is suitable also for the numerical implementation as 
it was illustrated on the particular example of the Ising 
dodecahedron. In addition, the method can be rather 
straightforwardly adapted in order to calculate magnetic 
properties of the irregular Ising spin clusters with a more 
complex geometric structure as well. 

It has been demonstrated that the magnetic behaviour 
of antiferromagnetic Ising spin clusters with the shape of 
regular polyhedra is basically influenced through a geo¬ 
metric spin frustration with exception of the Ising cube, 
which is the only bipartite regular polyhedron. In a conse¬ 
quence of that, the magnetization curve of the Ising cube 
involves just one trivial plateau at zero magnetization due 
to a simple antiferromagnetic ground state in contrast to 
the magnetization curves of all other regular Ising polyhe¬ 
dra. Among these, the Ising tetrahedron and icosahedron 
still exhibit a zero magnetization plateau, which is succes¬ 
sively followed by the other intermediate one-half plateau 
in the former case and respectively, the intermediate one- 
third and one-half plateaux in the latter case. On the other 
hand, there is a lack of zero magnetization plateau in the 
magnetization curves of the Ising octahedron and dodec¬ 
ahedron, which display only a single intermediate plateau 
at one-third and one-fifth of the saturation magnetization, 
respectively. 

The intermediate magnetization plateaux emerging in 
a low-temperature magnetization process owing to a ge¬ 
ometric spin frustration reveal highly degenerate ground 
states, the microscopic nature of which was clarified along 
with their respective degeneracies. We have found con¬ 
vincing evidence that the higher number of magnetization 
plateaux does not necessarily imply a quantum nature of 
higher-order magnetization plateaux contrary to conven¬ 
tional expectations [36, 37], but it may be just a simple 
consequence of geometric spin frustration. It has been 
also evidenced that the magnetization jumps accompany¬ 
ing a transition between different magnetization plateaux 
manifest themselves in a giant magnetocaloric effect dur- 
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ing the adiabatic demagnetization. From this point of 
view, the Ising octahedron and dodecahedron have turned 
out as the most prominent geometrically frustrated spin 
clusters that would allow an efficient cooling down to low¬ 
est achievable temperatures. 
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